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Abstract 

We study the interplay between nonlinearity in exciton transport and trap- 
ping due to a sink site for the dimer and the trimer with chain configuration 
by a numerical integration of the discrete nonlinear Schrodinger equation. 

Our results for the dimer show, that the formation of a self trapped 
state due to the nonlinear coupling increases the life time of the exciton 
substantially. Self trapping can be enhanced by the sink for short times, 
but for long times it disappears. In the trimer consisting of a subdimer 
extended by a sink site exists a transition between states localized on the 
two sites of the subdimer before for larger nonlinear coupling self trapping 
on one site of the subdimer is observed. 

For large trapping rates the fear of death effect leads to an increasing 
life time of the excitation on both, the dimer and the trimer. The sink 
site is then effectively decoupled. We explain this effect using an asymp- 
totic theory for strong trapping and demonstrate it by direct numerical 
computation. 
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1 Introduction 



The purpose of this paper is to study the interplay between a coherent transfer 
regime of an exciton and two processes leading to the loss of the linear character 
of the exciton transfer, namely trapping of the particle and transport nonlinear- 
ity. The system under study is a moving quasiparticle interacting strongly with 
polarization vibrations. At one of the sites of the configuration trapping with a 
prescribed rate 7 can occur. 

Trapping of quasiparticles constitutes an important phenomenon in many 
areas of physics. In photosynthesis, for instance, an exciton in a harvesting 
antenna transfers its energy to a so called reaction center, where it can be trapped. 
Electron transfer processes then follow. 

Much work has been done on the exciton transfer theories in the last two 
decades. The main effort has been directed to obtaining equations which de- 
scribe coupled coherent and incoherent motion of the excitation. The corre- 
sponding theories include Pauli Master Equations (PME), General Master Equa- 
tions (GME), Stochastic Liouville equations (SLE) and Continuous Time Ran- 
dom Walk (CTRW) [1-26]. Only recently [17-28] the problem of a rigorous 
description, leading to positive occupation probabilities, of the trap (modelled 
as a so called sink) in the coherent and near coherent regimes of the excitation 
transfer has been solved using the SLE, GME and CTRW methods. 

Using the SLE [17-20] we have shown for more extended systems, that in 
presence of a sink, the excitation behaves as it would like to avoid the site influ- 
enced by the sink. This effect has been observed also in the case of a semiinfinite 
linear chain [21] . We called it the "fear of death effect". In this case the decay 
of the whole occupation probability P(t) with time is slower for larger 7. 

In particular, in a linear trimer with a sink introduced at one end and for a very 
large trapping rate 7, the rest (dimer) could be taken as decoupled. The exciton 
moves from the very beginning approximately in the dimer avoiding the place 
influenced by the sink. Analogously, in the hexagonal model of the photosynthetic 
unit the reaction centre seems to be almost decoupled from the antenna system for 
a very large trapping rate 7. Still there is no doubt that after a long time the whole 
occupation probability P(t) of the exciton disappears due to the trapping term. In 
our recent paper [27] we investigated in some simple finite systems (dimer, trimer) 
the interplay between the trapping and the interaction with a bath described in 
the framework of the SLE with the formal Haken Strobl parametrization. 

In the last few years much attention has been paid to the problem of nonlinear 
interactions in coupled exciton- vibration systems [30-29]. A possible approach 
to it is the so called discrete nonlinear Schrodinger equation (DNLSE, sometimes 
also called discrete self trapping equation) which, though not yet rigorously jus- 
tified has the advantage of allowing for a simple analytical treatment. In this 
equation besides the transfer matrix element V a nonlinear coupling constant \ 
discribing the site energy lowering due to quasiparticle formation is contained. 
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One encounters the discrete nonlinear Schrodinger equation in the context of 
numerous phenomena from various fields of physics. 

Although the DNLSE was already mentioned in early papers on the polaron 
problem (see e.g. [52]) and much physics has been extracted through numerical 
analysis [29] exact solutions are not known in general. However, it has been 
found recently, that considerable insight can be gained into essential physics of 
the system through the exact analytical solutions when the system has only two 
sites. This case of the nonlinear dimer has been investigated during the last view 
years from many points of view. Kenkre with his collaborators [30-42] and Szocs 
and Banacky [43,44] have tackled the problem with the occupation probability 
difference as the most important quantity in mind. The exact time evolution 
for the occupation probability difference was given in terms of Jacobian elliptic 
functions for arbitrary initial conditions. Setting the occupation of the first site 
initially to one the transition to a self trapped state was obtained by Kenkre et al. 
[30] in a dimer with x ^ 4V. However, one sided oscillations in the occupation 
difference appear for some initial condition already from \ ^ 2V. This can be 
related to the appearence of a homoclinic structere on the Bloch sphere, which 
was shown by Esser and Hennig [45-49]. For increasing x the separatrix of this 
homoclinic structure grows and at \ = 4V starts to include the poles of the 
sphere which correspond to initially localized states. 

The limiting behaviour for t — > oo has been treated with and without dissipa- 
tion which was modelled in the framework of the stochastic Liouville equations 
with diagonal and nondiagonal Haken-Strobl parameters. Attention has been 
paid also to the influence of symmetry in the dimer, trimer and n-mer. 

Excitation transfer in trapless dimers and trimers can be described by Hamil- 
tonian flows located on the surface of the Bloch sphere for the dimer and its 
higher dimensional analog for the trimer case. This allows to use the methods 
of nonlinear dynamics in the analysis of the transfer properties, in particular to 
prove the appearence of bifurcations in the flow line picture [29,45-49]. Esser 
with his collaborators pointed out another interesting feature: the possibility of 
chaotic behaviour due to a perturbation of the homoclinic structure which con- 
stitutes a route to chaos. In this case one has to leave the description of the 
problem in terms of occupation probabilities only and must include the evolution 
of all Bloch variables, i.e. one must follow precisely the time development of the 
nondiagonal elements of the corresponding density matrix. In this connection 
interesting results have also been obtained in the nonlinear dimer, when the adi- 
abatic assumption is dropped, i.e. when vibrational relaxation which leads to the 
nonlinear transport is not assumed to be infinitely fast [50]. 

The aim of our contribution is to investigate the interplay between the nonlin- 
earity and the trapping of the excitation (here modeled as a sink). We shall follow 
the influence of both the nonlinearity and trapping on the transfer properties. 

The paper is organized as follows. In section II the formulation of the prob- 
lem and basic equations are given. Sections III and IV contain equations which 
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describe the dynamics of a dimer and a trimer, respectively, including nonlinear 
coupling and the influence of a sink placed at one end of the configuration. We 
present their numerical solutions and derive some analytic expressions for the 
case of very strong trapping. In section V conclusions are drawn. 



2 Formulation of the model and basic equa- 
tions 

2.1 The DST equation 

As the basic equation describing the excitation transfer in our systems we use 
the discrete nonlinear Schrodinger equation 

d 

^ ^ C-n (^) (^n X^l^l )c n ~\~ ^ ^ VnmC-m (1) 

where c n is the probability amplitude of the excitation (in what follows denoted 
as exciton) to occupy the molecule at the n-th site, e n the exciton site energy, 
Vnm {n 7^ rn) the transfer matrix element and Xn the nonlinearity parameter 
describing the lowering of the site energy by exciton occupation. Introducing the 
density matrix elements 

Pmn(t) ■= c m (t)c* n (t) (2) 
we describe the system evolution by the density matrix equation 
d 

«^|Onm (^n ^m) Pnm ~\~ ^ ^ (VnkPkm Pnk^km ) • (3) 

Here is 

^■n • &n }CnPnn (4) 

an effective site energy. 



2.2 Influence of the sink 

Trapping at a particular site s, consistent with the condition p nn (t) > for any 
time t, can be introduced into eq. (3) by supplementing the r.h.s. with the terms 

-i( 7 /2)p ns (n^s) (5) 

for the non-diagonal and 

-ilPss (6) 

for the diagonal matrix elements connected with the trap site s. 

The description of the influence of the trap as given by the terms (5), (6) is 
equivalent to introducing into the Hamiltonian the matrix element H ss = —i^f/2 
for the sink place s. 
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2.3 Perturbation theory for a large 7 

In the case of a large trapping rate 7 at site s, i.e. when 7 exceeds by far all the 
other energies entering the r.h.s of the eq. (3), one can simplify the system (3- 
6) by using quasistationarity in the equations for matrix elements p mn involving 
the trapping site. This corresponds to neglecting times smaller than 7 _1 in the 
evolution of the system and is formally achieved by setting p mn = for m = s 
or n = s (m ^ n) and m = n = s. This procedure is quite analogous to that 
used for the transition from the stochastic Liouville to the Pauli Master equation 
in [3,4] when quasistationarity is assumed for the non-diagonal density matrix 
elements due to the fast phase relaxation connected with the stochastic sources. 
Here we apply this approximation to a part of the nondiagonal density matrix 
elements, namely those connected with the trap site. As we will show, this results 
in analytical expressions for the fear of death effect in the case of large 7. One 
finds for the nondiagonal elements n ^ s 

Pns = ~, T7^y2(PnkVks — V n kPks) (7) 

t n - t s - z(7/2) ^ 
and for the diagonal elements connected with the trap 

Pss = - Y]{PskVks - VskPks)- (8) 

7 k 

From the r.h.s. of eq. (7) it is evident that the p ns are of the order of the 
small parameter (V/7) <C 1 where V is the typical value for the transfer matrix 
elements entering the r.h.s. of (7). Evaluating (7) and (8) in this small parameter 
iteratively one obtains approximations to p ns and p ss for the case of large 7 which 
complement our numerical results and are considerd below. We shall now treat 
the dimer and the trimer separately. 



3 The dimer with a sink 
3.1 Equations of motion 

Restricting the transfer to the two sites of a dimer with V i2 = V21 = —V < one 
arrives at 



ip\i 


= +V( Pl2 


- P21) 






IP22 


= -v( Pl2 


- P21) 


- 11P22 




iPl2 


= +V( Pll 


- P22) 


-1(7/2) 


P\2 + (ei - £2)^12 


iP21 


= -v( Pll 


- P22) 


-1(7/2) 


P2\ - (ei - £2)^21 
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Passing to the Bloch variables 



X\ := P12 + P21 

x 2 ■= i(pi2 ~ P21) (10) 
x 3 ■= P11 - P22 

and introducing the norm as a separate variable 

n := pn + P22 (11) 

one obtains the equations of motion 

xi = — (7/2) x 1 + Ae(t)x 2 

x 2 = -(7/2) x 2 - Ae(t) Xl + 2Vx 3 

x 3 = -( 7 /2)(x 3 -n)-2Vx 2 (12) 

h = -(7/2) (n - x 3 ) 

with Ae(t) := (e 2 — £1) + xx 3 . 

3.2 Previously obtained results 

In the absence of trapping, i.e. for 7 = 0, the system (12) reduces to the nonlinear 
Bloch equations considered in [45]. In this case besides the energy there is a 
second integral of the motion restricting the Bloch variable to the surface of a 
sphere. 

In [45] the fixed points of the trapless system (12) were analyzed. One finds 
that for x = 2V a bifurcation is realized which follows for the symmetric dimer 
(ei = £2) from x 2 = 0, i.e. 

(2V- X xi)x 3 = (13) 

For x ^ 2V the only solution of eq. (13) is given by x 3 = 0, whereas for x > 2V 
the new solution X\ = appears, resulting in two fixed points with x 2 = 0, 

x 3 = ±^1 — (~) 2 - These points correspond to the new ground states arising 
from the polaronic effect which holds the quasiparticle preferentially at one of 
the two sites of the dimer. The solution x 3 = for x > 2V corresponds to 
an unstable hyperbolic point at the center of a separatrix. Into this homoclinic 
structure the new ground states are embeded [45-47]. 

Kenkre et al. were able to rewrite the system (12) into a closed equation for 
the probability difference x 3 = pit) = Pi(t) — P 2 {t). The analytical and numerical 
solution leads for the case of the localized initial condition -Pi(O) = 1 to the so 
called self trapped states for x ^ 4V [30]. 

We now turn to investigating the influence of the nonlinearity on the trapping 
process of the exciton in the dimer. We shall check carefully changes in the 
dynamic properties due to the trapping rate 7. 
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3.3 Analytical results for strong trapping 

From (7), (8) we know that strong trapping leads to the expressions 

V 

Pl2 = 7-—-{p22-pXx) (14) 

e 1 - e 2 - i{j/2) 



[P21 = P12) 



an 



d 



iV 

P22 = [Pl2 ~ P2l) (15) 

7 

from which an explicit expression for the occupation of the sink site p 2 2 can be 
derived. In the case without nonlinearity (x = 0) one finds 

"° 2 = (e ,-J+(7/2)» (16> 

which demonstrates the fear of death effect: due to the strong trapping the 
occupation of the sink site is very small, namely p 2 2 ~ {V/l) 2 Pn *C Pn- The 
effect of nonlinearity in the perturbation series for strong trapping is to modify 
the energy difference in the denominator: 

V 2 

(ei - e 2 - XPn) 2 + (7/2) 2 
3.4 Numerical results 

We have calculated the time dependence of the Bloch variables in presence of non- 
linearity and trapping by a direct numerical integration of the nonlinear equation 
(I) including the trapping as described in Sec. 2.2. 

Our results are displayed in figs. 1-7 for different trapping rates 7 and 
various nonlinearity parameters x- We have chosen for all figures a symmetric 
configuration (e,\ = £2) with V = 1 and the initial condition -Pi(O) = 1. 

In fig. I the formation of a self trapped state for x ^ 4V is shown. The 
occupation difference oscillates between x 3 = — 1 and I with mean value zero for 
X < 4V. For x > 4 it remains in the vicinity of x 3 = 1, i. e. the exciton is 
preferentially located at the site I of the dimer. For exactly \ = 4V the initial 
condition -Pi(O) = 1 is located on the separatrix between trapped and detrapped 
states and approaches a very unstable hyperbolic fixed point [45] (due to the 
numerical inaccuracies the curve actually shown corresponds to x a little bit 
below 4.0). 

For both, 7 and x different from zero we show in the parts (a) of the figs. 2 - 
6 the time dependence of the total occupation probability nit) and in the parts 
(b) the occupation difference x 3 (t) rescaled with the total occupation probability. 

In fig. 2 the influence of the small trapping rate 7 = 0.5 is displayed. The 
total occupation probability decays slowly. Strong nonlinearity (x > 5) leads to 
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an effective "switching off" of the trapping. In this case the probability remains 
almost constant until t > 10 /V. 

In figs. 3 and 4 the medium trapping rates 7 = 1., 2., respectively, have been 
chosen. For the linear dimer \ = the life time of the exciton is smaller for 7 = 2 
than for large trapping rates 7 = 5, 10 (figs. 5, 6). This is due to the "fear of death 
effect" which becomes pronounced for 7 > 2. It leads to a very large difference 
in the occupation probabilities of the sites 1 and 2. In fig. 7 the influence of 
different trapping rates (indicated at the right) on the occupation probability at 
the trap site /^(O are presented for the linear dimer. /^(O becomes smaller for 
larger trapping rate and the decay with time of the whole occupation probability 
is consequently slower. 

A concise inspection of the probability difference P\(t) — ^2(^)5 sca led to the 
whole occupation probability (figs. 2-6 (b)) reveals, that for growing trapping 
strength the transition to the self trapped state takes place for \ < 4V already 
(but compare Fig. 2(b) in [32]). It is formed at the beginning of the time 
development and after some time it disappears. The self trapped state is not 
stable because the whole probability decays completely Pit) — > for t — > 00. 
Self trapping leads to a considerable enhancement of the life time of the exciton, 
e. g. for medium trapping rates (figs. 3, 4) and strong nonlinearity (x = 6) it is 
approximately four times longer than in the linear case x = 0. 

A very strong trapping rate destroys the influence of the nonlinearity almost 
completely (figs. 5, 6). The transfer probabilities for the exciton are reduced 
so much, that the occupation probability /^(O is very low and the decay of 
the whole occupation probability Pit) is almost independent of the nonlinearity 
parameter x- However, the formation of a self trapped state is still reflected in 
the normalized occupation probability for short times. 

4 The trimer with a sink 
4.1 Previously obtained results 

In the past some attention has been given to the solution of the nonlinear 
Schrodinger equation for the trimer. Preferentially symmetric trimers with sym- 
metric initial conditions have been studied. 

Eilbeck with his collaborators investigated the stationary solutions of the 
symmetric cyclic trimer. Their results are presented in Fig. 2 in [29]. They 
derived a bifurcation point and discussed the stability of the obtained solutions. 

Kenkre with his collaborators were able to obtain an equation for the proba- 
bility difference x 3 = pit) = P\(t) — 2P 2 (t) for a symmetric trimer with symmetric 
initial conditions [39,40]. 
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4.2 Equations of motion 

We consider a trimer with V\ 2 = V21 = — V and V23 = = — V < 0. A 
sink is introduced at site 3. Similiarly to the case of the dimer one can pass 
to generalized Bloch variables Xi_ 8 containing also the variables Xi_ 3 which were 
previously introduced for the dimer. The procedure of introducing the generalized 
variables Xi_ 8 is explained in the appendix. One obtains 



X\ 


— P12 + P21 


x 2 


= H/°12-/°2lj 


x 3 


= Pll - P22 


X4 


= P31 + Pl3 


x 5 


= i(pi3 - P3l) 


x 6 


= P23 + P32 


x 7 


= 1{P23 - P32) 


x 8 


= Pll - P33 


x g 


= P22 — P33 


n 


= Pll + P22 


N 


= Pll + P22 + P33 



Besides the linearly independent variables we have also introduced Xq — Xg — X3 as 
well as the occupation sum of the sites 1 and 2 n and the total trimer occupation 
N. The variable n will be useful for considering the sink influence at site 3. We 
obtain the generalized system of Bloch equations 

xi = +(e 2 — ei)x 2 — V'x 5 

x 2 = +(ei — e 2 )x 1 + V'x^ + 2yx 3 

x 3 = —2Vx 2 + V'xr 

= +(e 3 - ei)x 5 + Vx 7 - V'x 2 - (7/2) x 4 

x 5 = +(ei - £3)2:4 - Vx 6 + V'x t - (7/2) x 5 

x 6 = +(e 3 - e 2 )x 7 + Vx 5 - (~(/2) x 6 (19) 

X7 = +(£2 - £3)^6 - Vx 4 + 2V'x 9 - (7/2) x 7 

x 8 = —Vx 2 + V'x 7 — (^/3)(N — x 8 — x 9 ) 

x 9 = +Vx 2 + 2V'x 7 - (7/3) (N - x 8 - x 9 ) 

h = —V'x 7 

N = -( 7 /3)(N -x 8 -x 9 ) 



9 



4.3 Analytical results for strong trapping 

According to (7), (8) strong trapping leads to 



Pi3 = -(V> 23 - V' Pl2 ) (20) 

ei - e 3 - 27 

P23 = 1 V ,^ ( V 'P™ ~ V> P^ + V P™) ( 21 ) 

£2 - £3 - H7/^J 

and 

iV 

P33 = {P23 - P32)- (22) 

7 

This means that all density matrix elements with one index 3 are of the order 
0(V/"f) and therefore small, p 33 is even of order 0(V 2 /~{ 2 ). Keeping only terms 
of leading order one can express x 4 _7 through Xi_ 3 and obtains a closed system 
of equations for the latter: 

xi = -(7/2) xi + (e 2 - t 1 )x 2 

x 2 = -(- f /2)x 2 -(e 2 -€ 1 )x 1 + 2Vx 3 (23) 

x 3 = -(j/2)(x 3 -n)-2Vx 2 

h = -(7/2) (n - x 3 ) 

Comparing this to (12) we notice that the equations of motion for the trimer 
reduce for very strong trapping to the dimer equations with a newly defined sink 
rate 

== <ZZ. (24) 

7 

This eq. (24) once again demonstrates the fear of death effect: The effective sink 
rate 7 becomes small compared to V for 7 ^> V, i. e. the sink site in the trimer 
is effectively decoupled. 



4.4 Numerical results 

We have calculated numerically the time dependence of the site occupation prob- 
abilities in the presence of trapping and of a nonlinearity. With respect to the 
nonlinearity we consider two different cases: 

Configuration (I): The nonlinearity \ i s the same for all three sites 

Configuration (II): The sink site is linear, i. e. Xi = X2 = X? X3 = 

Our data were obtained through a numerical integration of the nonlinear Schro- 
dinger equation (1). We have chosen for all figures V = 1 and the initial condition 
-Pi(O) = 1. The sink influences site 3. In the analysis we will concentrate on self 
trapping, the fear of death effect and regular vs. chaotic time dependence. 
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For reference we show first the linear trimer without trap fig. 8a. In this case 
the exciton coherently oscillates between the ends of the trimer. The sites 1 and 
3 have equal average occupation probabilities which are due to the special initial 
condition -Pi(O) = 1 larger than the occupation of site 2. The situation is similiar 
in the case of a small trapping rate 7 = 0.2 (fig. 8b) but the total occupation 
probability in this case slowly decays. Increasing the trap strength up to an 
intermediate value 7 = 2 (fig. 8c) makes the decay faster until for 7 larger than 
~ 2 the fear of death effect becomes pronounced and the life time of the exciton 
grows again but now with the occupation probability at site 3 much smaller than 
for small trapping rates. We show for example the situation at 7 = 10 (fig. 8d) 
when the exciton oscillates between the sites 1 and 2 and the trap site is almost 
decoupled. 

A moderate nonlinearity (up to x ~ 2) does not change the appearance of the 
graphs very much (fig. 9). Again the fear of death effect controls the behaviour 
for larger trapping rate and in this case there is not much difference between the 
two configurations (I) and (II, not displayed) because the occupation of site 3 is 
reduced so much. In contrast, for weak trapping the regular oscillations of type 

(I) change into a quite irregular behaviour as one removes the nonlinearity at site 
3 (figs. 9a and b). This effect is also present without the trap (not displayed) 
and is obviously due to the destruction of the symmetry between the sites 1 and 
3. 

As one increases the nonlinearity further, a significant change of the behaviour 
towards a self trapped solution occurs. We discuss this first for the symmetric 
trimer of configuration (I) without trap. The regular behaviour up to x = 2 (fig. 
10a) is replaced by a rather chaotic one for x = 3 (fig. 10b). Still the average 
occupation probabilities are approximately the same for all three sites. There 
is, however, a large intervall where the occupation of site 3 is very small. For 
X = 4 (10c) the transition to a self trapped state on the dimer has occured and 
the occupation of site 3 remains small all the time. The appearance of the time 
dependence is much more regular than close to the transition point in fig. 10b. A 
further increase of the nonlinearity parameter causes the exciton to be completely 
localized at site 1 (fig. lOd) for x ^ 5, after that no other qualitative changes 
occur, just frequency and amplitude of the small oscillations in the occupation 
probabilities slightly depend on x ( n °t displayed). The consecutive transitions 
can be explained by the bifurcations which the stationary solutions undergo as the 
nonlinearity changes. They have been discussed in [29] in detail for the nonlinear 
Schrodinger equation without sinks. 

The following figure 11 shows that the bifurcations are still present for a small 
trap strength 7 = 0.2 (cmp. also figs. 8b; 9a, b). In the symmetric conf. (I) the 
transition to self trapping on the dimer occurs for all sink rates later than in conf. 

(II) . We illustrate this with the figs. 11a, b. The effect can be understood having 
in mind that equal site occupation probabilities for the asymmetric configuration 
(II) cause a shift in the effective site energies e 8 - on the dimer which in turn 
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obstructs the transition from the dimer to site 3. The transition to self trapping 
on site 1 is not affected by the two different configurations (not displayed). 

Even for very large 7 there is a clear transition from self trapping on the 
dimer to a solution completely localized at site 1 (figs. 8d; 9d and 12). In this 
case, however, the occupation probability on site 3 is always small. This is now 
not due to nonlinearity induced self trapping but to the fear of death effect. In 
general, both effects work together in slowing down the decay of the excitation. 
The life time is e. g. for 7 = 10 and x = 5 doubled in comparison to x = 3. 

The formation of self trapped states due to nonlinearity is not stable for 7 > 0: 
as the the total occupation probability becomes small with increasing time the 
effect disappears and the site occupation probabilities level out, as it can be seen 
e. g. from fig. lib it > 35) and fig. 12b it > 30). The nonlinearity effects must 
decrease for large time since the occupation probabilities entering the effective 
site energies (4) then become negligibly small. 

In contrast, the lowering of the occupation p 33 due to strong trapping does 
not depend on time. The sink site is almost decoupled and the behaviour of Pi 
and P 2 can be interpreted as a perturbed dimer (figs. 9d and 12a, b). We note the 
interesting effect, that the time dependence looks quite regular as well for X ^ 3 
(fig. 9d) as for X ^ 5 (fig. 12b) but much more chaotic for x = 4. We relate this 
to the fact, that for x = 4 the initial condition -Pi(O) = 1 is just on the homoclinic 
orbit (cmp. fig. 1) which encloses the self trapped states for the system without 
trap. Time evolution will therefore take the system in the perturbed dimer very 
close to the hyperbolic fixed point where it is very unstable and as a result one 
observes chaos (see [45] - [49]). 

In the last fig. 13 we present a numerical check for eq. (24). We compare 
the time evolution for a trimer with the large sink rate 7 = 10 (full lines) to a 
dimer with 7 = 0.4 (dots) according to this equation. We find good agreement 
in particular for x = (a). For x = 3 the result is still reasonable, but not as 
good as in (a) since (24) was derived also under the assumption 7 ^> x- 

5 Conclusions 

In the present paper we studied the combined effect of trapping due to a sink 
and transport nonlinearity due to excitonic - vibronic interactions on the transfer 
dynamics and the life time of an exciton. The specific systems we have analyzed 
were the dimer and the trimer with a chain topology. A sink was introduced 
at one end of the configuration using the formalism given in [28]. The initial 
condition for the numerical integration of the NLSE was always a state localized 
at the end opposite to the sink site. 

The main effect of the nonlinearity is self trapping. In the dimer with sink 
a tendency to form a self trapped state could be observed for x < 4V, i.e. for 
smaller nonlinearity than in the isolated dimer. In the trimer self trapping on 
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the sinkless subdimer for an intermediate nonlinearity was found besides states 
which are localized on one site only for strong nonlinearity. The transition to 
a self trapped state was accompanied by relatively irregular time dependencies. 
Though self trapping was enhanced by the sink for short times the decay of the 
total occupation probability destroyed the localized states for long times. Self 
trapping was shown to have a crucial influence on the life time of the excitation 
which grows for increasing nonlinearity. 

For large sink rates the fear of death effect controls the decay of the excitaton. 
In this case the occupation of the sink site is reduced so much that the life time 
even increases with growing sink rate. The sink site is then almost decoupled 
from the rest as it was shown using an asymptotic approximation. 

For small trapping rates there was a clear change in the behaviour when the 
nonlinearity of the sink site was was dropped while this effect was supressed by 
the effective decoupling of this site for strong trapping. 
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Appendix: Generalized Bloch representation for 
the trimer 



The Bloch variables for the dimer are derived from the state vector c 
by employing the Pauli matrices Oi 



Cl, C 2 



1 

1 



0"2 



-i 

1 



0"3 



1 
-1 



(25) 



which are the generators of the Lie algebra SU(2). In this case one defines 

(i = 1, 2, 3). (26) 



C + OiC 



Analogously, one derives generalized Bloch variables for the trimer from the state 
vector c = (ci,c 2 ,c 3 ) T using the Gell - Mann matrices A;, which are generators 
of the Lie algebra SU(3), in the form 



Ai 
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They generate the variables 



c + A,C 



1, 



9). 




(27) 



(28) 



These variables X{ are by definition real. We have chosen A 8 in a modified form 
convenient for our calculations and added A 9 which is linearly dependent (A 9 = 
A 8 — A 3 ). Our variables x 8 and x 9 are therfore generalizations of the variable x 3 
and have the form of differences between occupation probabilities. 
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Figure 1: Occupation probability difference x^(t) for the trapless dimer (j = 0) and 
various nonlinear coupling strengths x (given to the right of the plots). The transition 
to a self trapped state at site 1 (x% ~ 1) is obvious. 
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Figure 2: Total occupation probability n(t) (a) and normalized occupation probability 
difference x^(t)/n(t) (b) for a dimer with weak trapping 7 = 0.5 and various nonlinear 
coupling strengths x (given to the right of the plots). Again the formation of a self 
trapped state for x > 4 can be observed. It is however not stable and decays after some 
time (t ~ 10 for x = ^) 




Figure 3:Total occupation probability (a) and normalized occupation probability dif- 
ference (b) for a dimer with intermediate trapping 7 = 1.0. The decay of the total 
occupation is faster than for 7 = 0.5 and also the self trapped state disappears earlier. 
The influence of self trapping on the life time of the exciton is evident. 
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Figure 4: Total occupation probability (a) and normalized occupation probability differ- 
ence (b) for a dimer with intermediate trapping 7 = 2.0. The tendency to form a self 
trapped state can be recognized even for x = 3. 
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Figure 5: Total occupation probability (a) and normalized occupation probability dif- 
ference (b) for a dimer with strong trapping 7 = 5.0. The life time of the exciton is 
similiar to the case 7 = 2 but clearly shorter than for 7 = 10. The oscillations of the 
occupation difference have disapeared and the occupation at site 1 is always much larger 
than at the sink site 2 (fear of death effect). 
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Figure 6:Total occupation probability (a) and normalized occupation probability differ- 
ence (b) for a dimer with very strong trapping 7 = 10. The exciton is almost completely 
localized at site 1. The influence of the nonlinear coupling is very small. 
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Figure 7 '.Occupation probability P 2 (t) of the site with the trap for no nonlinearity 
(x = and various trapping strengths 7 (given to the right of the plots). The 
fear of death effect leads to the very small occupation of the sink site for 7 > 5. 
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Figure 8:Decay of an exciton on a linear trimer (x = 0) for different 
sink rates: (a) 7 = 0, (b) 7 = 0.2, (c) 7 = 2, (d) 7 = 10. The fear of 
death effect causes the life time of the exciton to be shorter for 7 = 2 
than for 7 = 10. 
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Figure 9:Decay of an exciton on a trimer with small nonlinear coupling 
X = 2.0 for different sink rates and configurations: (a) 7 = 0.2, conf. 
(I) X = Xi(i = 1, A 3); (b) 7 = 0.2, conf. (II) X i = X2 = X, Xs = 0; (c) 
7 = 2, (d) 7 = 10, both conf. (I). Conf. (II) shows irregular behaviour 
for small sink rate due to the destruction^ symmetry (b). For larger- 
sink rate conf. (II) is not displayed because it behaves very similiar to 
conf. (I). Comparison of (c) and (d) again shows the fear of death effect. 



(a) (b) 




5 10 15 20 25 30 35 40 45 50 5 10 15 20 25 30 35 40 45 50 
t t 

(c) (d) 



\/\ 


























a. 


/\A/\AA/V\A/V\ArV\A/^A/W\/V\A/v\ 























5 10 15 20 25 30 35 40 45 50 5 10 15 20 25 30 35 40 45 50 
t t 



Figure 10:The formation of a self trapped state for the trapless sym- 
metric trimer (conf. (I), ~f = §): At x = 2 (a) the behaviour is very 
similiar to the linear trimer of fig. 8a. The time dependence becomes 
irregular for x = 3 (b). For x = 4 (c) the solution is self trapped on a 
dimer and for x = 5 (d) on site 1. 
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Figure ll:Formation of a self trapped state for a trimer with the small 
sink rate 7 = 0.2: (a) x = 3, conf. (I), (b) x = 3, conf. (II), (c) x = 4, 
and (d) x = 5, both conf. (I). For conf. (II) self trapping on the dimer 
occurs for smaller nonlinearities. The self gapped state on the dimer in 
(c) dis s appears for t > 35. 
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Figure 12: Transition to a self trapped solution on site 1 for strong 
trapping 7 = 10. The nonlinear coupling is (a) x = 4 and (b) x = 5. 
The self trapped state is unstable for / > 30. The occupation of site 3 is 
always very small (fear of death effect). The time dependence for x = 4 
is very irregular. 
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Figure 13: Reduction of the trimer with strong trapping to an effective 
dimer with weak trapping. The full line shows the occupation probabilities 
of the sites 1 and 2 for a trimer with 7 = 10 and (a) x = 0, (b) x = 3. 
The dots show according to eq. (24) the corresponding effective dimer 
with the sink rate 7 = 0.4 
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